Study of spin-Peierls transition in a'-NaV^Os by infrared reflectivity 
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Polarized infrared reflectivity measurements have been performed on single crystals of the spin- 
Peierls compound a'-NaV20s in the temperature range 20-300 K. Pronounced spectral features 
associated with the formation of the dimerized phase were detected both in the a- and b-polarizations 
(perpendicular and parallel to the spin-1/2 chains, respectively). The temperature dependence of 
a salient spectral line at 718 cm -1 sharply rising below the transition temperature Tsp obeys a 
(1 — T/Tgp) 2 ' 3 law with T$p — 34.3 K and (3 ~ 0.25. In addition, a continuum signal is observed in 
the whole temperature range in the a-polarized optical conductivity spectra. In order to interpret 
these results, calculations of the static dimerization and of the optical conductivity based on a 
mean-field and a dynamical treatment of the lattice respectively are proposed. 
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The spin-Peierls (SP) transition has become a re- 
vived challenging effect in quasi-one dimensional mag- 
netism since the recent discovery of two inorganic tran- 
sition metal oxides compounds, first copper germanate 
CuGe0 3 (Ref. @) and more recently sodium vanadate 
a'— NaV205 (Ref. ||). These compounds contain linear 
spin-1/2 chains coupled to three-dimensional phonons. 
At low temperatures this coupling leads to the SP phase 
transition characterized by a lattice dimerization to- 
gether with the opening of a spin gap. 

The revival of the subject already investigated before 
in numbers of organic salts™ is due to the availability of 
large size and high quality crystals offering the-possibility 
to carry out a number of relevant experimentsu including 
x-rays and neutron scatteringcTEl. 

Optical spectroscopy is known to be a sensitive, tool 
to investigate properties of magnetic compoundsB and 
provides a detailed knowledge of phonons, magnons and 
their fluctuation states. IR transmission and reflectance 
experiments have been used to study the undistorted and 
SP phases in CuGe0 3 (Ref. §^U]). The SP phase of 
NaV2C>5 has been recently studied by IR absorption in 
the frequency region 50-400 cm -1 . New phonon-lines 
have been found to accompany the SP transitiont2lll3. 
Recently, the activation of a zone boundary phonon in the 
SP phase of CuGe03 has been observed as a new single 
tiny line in low temperature reflectivity spectraHI. A dis- 
tinctive feature of NaA^Os in comparison with CuGeC>3 
is the negligible magnetic interaction between the chains 
because of their isolation by V 5+ 05 non magnetic chains. 
Frustration ratios a w 0.24-^0.36 (Ref. Q|5|) have been 
suggested for CuGe0 3 and a m for NaV 2 5 (Ref. |). In 
addition, the nearest-neighbor intrachain exchange con- 
stant J~440 K, the transition temperature Tsp= 35 K 



and the spin gap A sp j n = 79 cm -1 in NaV^Os are typi- 
cally 3-4 times larger than in CuGeOa (Ref. |l3|). Also the 
dimerization 5, estimated in the adiabatic approximation 
from the T = spin gap, has the value 6 ~ 0.014 (Ref. 
0) and 8 ~ 0.048 (Ref. fjj|) for CuGe0 3 and NaV 2 5 , 
respectively. Thus, the magneto-elastic coupling plays a 
larger role in NaV20s and enhanced spectral signatures 
of this coupling are expected in this system. 

In the following, we report on infrared reflectivity mea- 
surements of NaV205 single crystals over frequencies be- 
tween 100 and 10000 cm" 1 and for temperatures between 
20 and 300 K. The study focuses first on the room tem- 
perature display of the excitation spectrum for the inci- 
dent polarization along the a- and b-axes (perpendicular 
and parallel to the spin-1/2 chains, respectively). Then, 
prominent new lines giving salient evidence of the dimer- 
ized phase below 34 K are reported. A detailed study 
of the temperature dependence of the line at 718 cm -1 
is presented. Finally, the outline of a calculation using 
Exact Diagonalization (ED) and Quantum Monte-Carlo 
(QMC) algorithms based on a mean-field or a dynamical 
treatment of the lattice is given. 

High-quality single crystals NaV205 with size up to 
2.7x5.1x0.3 mm 3 were prepared according to the pro- 
cedure described in Ref. [l7| The optical measurements 
were done with a Bruker IFS 113V spectrometer at nearly 
normal incidence of infrared radiation polarized along the 
a- and b-axes and k || c. Room temperature transmis- 
sion and reflection infrared spectra were measured in a 
spectral range between 90 and 10000 cm -1 with a reso- 
lution ranging from 1 cm -1 in the far- infrared region to 
5 cm -1 in the near-infrared region. The low temperature 
far-infrared spectra were measured with a continuous He 
flow cryostat down to 20 K with the absolute accuracy 
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FIG. 1. Room temperature reflectivity spectra and opti- 
cal conductivity of NaV^Os single crystal for a-axis (a) and 
b-axis (b) polarizations. 
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FIG. 2. Detailed temperature dependence in the reflectiv- 
ity of a'— NaV20s single crystal. (a)-E || a, temperatures 
from up to down: 25, 29, 32, 37, 300K. (b)-E || b, temper- 
atures from up to down: 21.8, 24, 27.1, 29.5, 32, 33.5, 34.3, 
35.3, 38, 300K. Spectra are vertically shifted by 7% (a) and 
10% (b). 



of temperature control about 0.1 K. The spectral range 
was 90-700 cm -1 for a-polarization and 90-1000 cm -1 
for b-polarization with the resolution 1 cm -1 . 

We can distinguish three different regions in the mea- 
sured infrared spectra of NaV2C>5 at 300 K. A rich 
phonon structure at frequencies below 1000 cm -1 is fol- 
lowed by transmission windows at 1500-5150 cm -1 for 
a-polarization and at 1120-5320 cm" 1 for b-polarization. 
At higher frequencies, a wide maximum centered at ap- 
proximately 7100 cm -1 was observed in the a-polarized 
reflectivity spectra. We have here restricted the discus- 
sion to the phonon part of the spectra. 

In Fig. [I] the room temperature a- and b-axes re- 
flectivity spectra (R) and the frequency-dependent con- 
ductivity (a) derived via Kramers-Kronig analysis are 
shown in the region of lattice vibration frequencies up 
to 1000 cm -1 . Three phonons are detected for the a- 
polarization (loto ~ 524, 251, 145 cm -1 ) and four for 
the b-polarization (ujto ~ 581, 367, 229, 175 cm -1 ), fa- 
axis polarized reflectivity spectra have a typical form of 
phonon bands in an insulator which manifests in dynam- 
ical conductivity as Lorentzian peaks. Obvious features 
in the a-polarized a spectrum are the occurrence of a 
broad excitation band ending at 400 cm -1 and an asym- 
metric line at 524 cm -1 with an underlying continuum 
extended above 1000 cm -1 (Fig. |l]a). However, since the 
line shape depends strongly on the high-w approximation 
of the data (> 10000 cm" 1 ), reflectivity measurements in 
the visible and ultraviolet regions are necessary for a def- 
inite analysis. 



Below Tsp, a pronounced modification of the reflectiv- 
ity spectra is observed with sharply rising new features at 
363, 200 cm" 1 (E || a), 718 cm" 1 (E || b) and 138 cm" 1 
(the polarization is unknown because of a poor signal- 
to- noise ratio in the very far-infrared region) . The a-axis 
polarized broad band at oj < 400 cm" 1 is enhanced at low 
temperatures, but neither the position nor the intensity 
changes at T = T SP (Fig. |). 

We performed a careful investigation of the 718 cm" 1 
line in the conductivity over temperatures between 20 
and 38 K since this line is clearly associated to the lat- 
tice transformation at T < T$p (Fig. ||). At low temper- 
atures (T=21.8, 24, 27 K) this peak is well described as 
a single Lorentz line with a constant linewidth. On the 
other hand, above 29 K the line profile is well described as 
superposition of a temperature-dependent Lorentz peak 
on top of a temperature-independent asymmetric broad 
line which can be clearly seen at 35.3 K and 38 K. This 
broad line might well be the signature of structural fluc- 
tuations growing in the vicinity of T$p. Indeed, after 
subtraction of the temperature-independent part, a per- 
fect Lorentz line with constant linewidth (8.2±0.5) cm" 1 
is obtained for all T below the transition. We believe 
that the Lorentz-type line corresponds to a folded zone 
boundary phonon mode activated by the structural dis- 
tortion and, therefore, the intensity of this line should 
be proportional to the squared dimerization. The tem- 
perature dependence of the normalized line intensity is 
shown in Fig. as well as the result of a fit with the 
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FIG. 3. Temperature evolution of the peak at 718 cm 
(temperatures from up to down: 21.8, 24, 27.1, 29.5, 32, 33.5, 
34.3, 35.3, 38 K) observed for a b-axis polarization. Spectra 
are vertically shifted by 0.9 Q~ cm - . The insert illustrates 
the correction procedure, as described in the text (□: 34. 3K, 
: 38K, V: difference, line: Lorentzian fit). 



(1 -T/T S p) 2/3 law. This fit in the temperature range 30- 
35 K gives T S p=(34.3±0.2) K, /3=0.25±0.1. One should 
note that the relatively small error on /3 achieved with 
only few points near the Tsp is due to the correction pro- 
cedure described above together with the low dispersion 
of experimental data as evident in Fig. || (insert). 

We proceed now to discuss the above results in the con- 
text of existing theoretical approaches. The crudest de- 
scription of the spin-Peierls ground state (GS) is based on 
the dimerized Heisenberg model: the antiferromagnetic 
coupling between nearest neighbor spins is assumed to 
be statically modulated (see Eq. (^) below) along each 
spin-1/2 chain. In fact, the physical origin of such a mod- 
ulation relics on the dynamical magneto-clastic coupling. 
One of the supplest model including dynamical effects of 
phonons readli3L3 (omitting the chain index), 

H = J]T(1 + A^.S.+i + H % + i/ c i astic , (1) 

i 

where A is the magneto-elastic coupling constant and Ui 
correspond to atomic displacements which modulate the 
exchange integral J. Here, we restrict ourselves to a sin- 
gle polarization of the phonon modes. Assuming inde- 
pendent optical phonons, the phonon hamiltonian takes 

2 

the usual form, i?° h = J^iijti + ^ K \\ u i}- An inter- 
chain elastic coupling H^ astic is required to account for 
the three-dimensional coherence of the phonons and in- 
volves an additional elastic coupling K± between neigh- 
boring chains. H~\ mtic is needed to obtain a finite value 
of Tsp- Note that Kn and A depend strongly on the ac- 
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FIG. 4. Temperature dependence of the normalized 
squared dimerization obtained by ED (14 (0), 16 (o) sites), 
QMC (□, 40 sites) and intensity for the phonon mode ob- 
served at 718 cm -1 in b-axis polarization (•). The solid line 
is a fit to a (1 - T/T S p) 213 law. All data are normalized to 
their 22 K value. Error bars are shown for the QMC and 
experimental data. 



tual direction in space of the displacements Ui which is 
not explicitely specified in (|]) (Ref. pl| ). 

The optical conductivity is theoretically defined as the 
GS correlation a{ui) cx (A\ R (u) + iQ + — H)~ 1 Aip) where 
the operator Ajp — ^ . Ai describes the coupling be- 
tween the (q = 0) photons and the medium. Quite gen- 
erally, one expects Ai cx E a u i T a p(i) 1 where E is the 
electric field of the light, is the polarization of the 
phonon mode and T a p{i) is the IR tensor which might 
include local spin fluctuations^! 

In lowest order, Ajp corresponds to the emission of 
a q = phonon i.e. T a p(i) — 8 a p. According to Ref. 
[ij, at T > Tsp there are 15 E || a and 7 E || b phonons 
which should be active in the first order IR processes. We 
did not observe all the predicted modes probably because 
some of them have small oscillator strengths. The group 
theoretical analysis in the SP phase is not available yet. 

Calculated optical conductivity spectra (T = 0) ob- 
tained by ED of a N=12 site ring are shown in Fig. || 
for a typical phonon frequency Q — J and several val- 
ues of the dimensionless coupling constant g = \/\/2MQ, 
(-^clastic i s neglected here and a variational basis is used 
for the phononscJ). As shown in the insert, the spin gap 
(and hence the dimerization) increases strongly with the 
coupling g. The doubling of the unit cell leads to an 
additional peak in a(oj) corresponding to a new q = 
phonon mode at higher energy resulting from a folding 
of the Brillouin zone. This mechanism can then explain 
the experimental observation of the new line at 718 cm -1 
if one assumes that the atomic displacements ui have a 
component along the chains .direction (b-axis) . 

As-suggested theoreticallyE3 and observed experimen- 
tallyEil, higher order processes involving one phonon and 
two-magnon scattering can also take place with, e.g., an 
IR tensor like T a p(i) = Sj-S i+ i(<5 Q ^ — e^ l+1 eg l+1 ), where 
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FIG. 5. Optical conductivity vs frequency calculated 
within model (|l|) for Q/J = 1 and for various spin-phonon 
coupling strengths g. Insert: A sp j n vs g. 



e vi+i j s the urn t vector along the chain direction. In 
this case, the maximum absorption would occur for E 
polarized along a. A continuous background is indeed 
observed for such a geometry (see Fig. |l|(a)) and could 
be associated to a two magnon phonon assisted absorp- 
tion. More detailed theoretical and experimental inves- 
tigations of this effect are left for a future study. 

Although the treatment of the phonon dynamics at fi- 
nite T is beyond the scope of this work, a simpler mean 
field treatment of the lattice is justified if the relevant 
phonons have a pronounced 3D character (i.e. K± is not 
too small compared to Kn ) and if Q,/ J < 1. In the adia- 
batic approximation (which corresponds, strictly r speak- 
ing, to the limit fi/J — > 0) the hamiltonian readso, 
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MF 



1 



J}^(l + (-iyS)S i -S i+1 + -NK5" 



(2) 



where N is the number of sites, K an effective elastic con- 
stant and S = S(T) the static dimerization which depends 
on the system size. The value of S(T) is obtained by min- 
imizing the total free energy Ft — F spin (T, S) + ^NKS 2 , 
where F sp i n corresponds to the spin part of hamiltonian 
(J|). The zero temperature value (5(0) ~ 0.048 is im- 
posed by the actual valua of the ratio A spin /J ~ 0.2 
known from experimental!]. The value of 6(T = 0), in 
turn, enforces the value of the ratio K/J ~ 3. Note 
that K and J are the only two free parameters of the 
model (J ~ 440 K is-determined from a fit of the high- 
T spin susceptibilitytj) and, hence, the theoretical value 
for Tsp comes as an output of the model. The QMC cal- 
culation suggests Tsp ~ 33 K, a value very close to the 
experimental one. The explicit temperature-dependence 
S 2 (T) obtained by ED and QMC techniques^ agrees rea- 
sonably with the experimental data, as shown, in Fig. ^, 
apart from the value of the critical exponents. 

In summary, we have studied the infrared reflectiv- 
ity on NaV2C"5 single crystals at temperatures down to 
20K. The SP transition order parameter was measured 



through the phonon spectrum temperature dependence. 
These results are compared with calculations based on 
mean field or dynamical treatments of the lattice. 

We acknowledge A. Damascelli and D. van der Marel 
for communicating their results of a similar study prior 
to publication. 
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